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Quantum fc-core percolation is the study ol quantum transport on fc-core percolation clusters 

where each occupied bond must have at least k occupied neighboring bonds. As the bond occupation 

probabihty, p, is increased from zero to unity, the system undergoes a transition from an insulating 

. . . phase to a metallic phase. When the lengthscale for the disorder. Id, is much greater than the 

04 ' coherence length, Ic, earlier analytical calculations of quantum conduction on the Bethe lattice 

demonstrate that for k — 3 the metal-insulator transition (MIT) is discontinuous, suggesting a new 

^^ . universality class of disorder-driven quantum MITs. Here, we numerically compute the level spacing 

Cn ' distribution as a function of bond occupation probability p and system size on a Bethe-like lattice. 

^ ' The level spacing analysis suggests that for k = 0, pq, the quantum percolation critical probability, is 

Cy , greater than pc, the geometrical percolation critical probability, and the transition is continuous. In 

contrast, for k — 3, pq — pc and the transition is discontinuous such that these numerical findings are 

consistent with our previous work to reiterate a new universality class of disorder-driven quantum 

O ! MITs. 
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I. INTRODUCTION 



With the exception of transition metal compounds, there exist two conventional paradigms for metal-insulator 
transitions (MITs): a Mott-Hubbard-type transition and an Anderson- type transition. The former is a consequence 
of tuning the interactions between electrons by changing the distance between atoms, for instance, and the MIT is 
discontinuous [IJ. The latter is a consequence of the tuning the disorder in the material, and the MIT is continuous [3, 
y]. Recently, a discontinuous disorder-driven MIT has been predicted for a model in which there exists geometric 
constraints on the disorder [3]. The model has been dubbed quantum fc-core percolation. Quantum k = 0-core 
i-rt percolation, or quantum percolation (QP), has been studied since the 1980's and exhibits a continuous, Anderson- type 
(~| : MIT ^5, 6] . Quantum percolation is defined accordingly. Consider a lattice whose bonds are occupied independently 
O ' and at random with bond occupation probability p. An electron can only hop between lattice sites i and j along an 
O, I occupied bond and cannot hop along an unoccupied bond. In addition, there exists a constant on-site binding energy 
which is set to zero for convenience. In total, the tight-binding Hamiltonian for quantum percolation is 

>'. H = Y,Uj4a,+h.c. (1) 

G^ ■ 
^. 

'£j ■ in which, 



»j 



CO . + _ / -'^ with probability p , , 

O ■ ^-^ ~ 1 with probability l-p ^ > 



and a\ and a^ are electron creation and annihilation operators. 

Quantum percolation exhibits a MIT as p is increased from zero, at least for three dimensions and above [7H9|. 
In two dimensions, some studies indicate a transition |10l - [l4| . while others do not [1^, [l3|- Analytical work on the 
Bethe lattice indicates that the quantum percolation transition may be in the same universality class as geometric 
C^ I percolation except with the transition probability Pq above which there exists extended zero-energy wavefunctions 
is greater than the threshold above which there exists a spanning cluster, pc [l3,[l3- The geometric percolation 
transition is a continuous and, therefore, so is the quantum percolation transition presumably. 

Recently, we have studied quantum fc-core percolation on the Bethe lattice [1]. The term fc-core refers to a geo- 
metrical constraint where every occupied bond must have at least fc — 1 occupied neighboring bonds [l9l - l24| . Such 
a geometric constraint may have implications for glassy systems [25], jamming systems [23], and even biological sys- 
tems [2q |. To enforce this constraint, bonds are initially occupied independently and at random with probability p. 
Then, those occupied bonds with less than fc — 1 occupied neighboring bonds are rendered unoccupied. This removal 
procedure proceeds recursively throughout the lattice until all occupied bonds satisfy the fc-core constraint. See Fig. 
1 for an example on the Bethe lattice. For fc > 3, the geometric percolation transition is in a different universality 
from ordinary geometric percolation with, for example, the fraction of occupied bonds in the spanning fc-core cluster, 
Poo, scaling with p as Poo = Pq + Pi{p — PcY^"^ for p > Pc (with Pc — p « 1) where Pq and Pi are constants and 
Pc is the critical occupation probability for classical (geometric) percolation. This result is to be contrasted with 




FIG. 1: Here fc = 3 and z — 4. The shaded circles denote branches that are fc — 1 connected to infinity. The removal of bonds 
1 and 2 eventually triggers removal of bond 3 and bonds emanating from vertex a, including the shaded circles. The remaining 
three branches emanating from the center site survive the removal process. 



Poo oc (p — Pc) for p > Pc just above the transition for fc = 0, 1. In other words, the geometric percolation transition 
is a random first-order transition, i.e. a discontinuous transition with several diverging lengthscales. 

Once the fc-core constraint has been implemented, we then impose the quantum percolation property that electrons 
can only hop along occupied bonds and compute the quantum conduction through the system self-consistently after 
assuming that the electronic wavefunction randomizes between bonds, i.e. the lengthscale of the disorder. Id, is much 
greater than the coherence length, Ic- Within this scheme, for fc = 0,1, the model reduces to ordinary quantum 
percolation. In this case, we found that 1 > Pq > Pc inhere Pc signals the onset of the geometric percolation transition 
and Pq signals the onset of quantum conduction. We also found a random first-order transition for fc > 2, and 
that interestingly enough, the critical threshold is the same as fc-core geometrical percolation critical point. This 
transition should be contrasted with the Anderson-type MIT, which yields a continuous transition as does the quantum 
percolation transition on the Bethe lattice. 

In this work, we numerically investigate the level statistics of quantum fc-core percolation to compare with our 
previous analytical results on the Bethe lattice obtained in the limit Id » Ic- In other words, how robust are our 
previous results in identifying a new universality class of disorder-driven MITs? Level statistics and its roots in 
random matrix theory is an important tool for studying universality, for example [27| . Correlations between energy 
eigenvalues of an individual quantum particle in a random potential in the conductive regime agrees with results from 
Gaussian matrix ensembles [28i]. In the localized regime, correlations between the energy eigenvalues are absent and 
the level statistics become Poissonian. Right at the MIT, however, the level statistics are distinct from Gaussian 
matrix ensembles [231. While these results pertain to the Anderson model where the disorder is on-site, the same 
analysis has been applied to quantum percolation on a cubic lattice, where the disorder is off-diagonal, and similar 
results have been found [3y|. In fact, the critical exponent for the divergence of the localization length extracted from 
the level statistics analysis is somewhat consistent with the Anderson model. Therefore, we will implement a similar 
numerical analysis for quantum fc-core percolation on a Bethe- like lattice to go beyond our previous approximation [J] . 

We, however, will not investigate the level statistics of quantum fc-core percolation on finite-dimensional lattices for 
now. It turns out that fc-core percolation on finite-dimensional lattices either exhibits properties of fc = geometric 
percolation or no transition. See, for example, Ref. [31] . One has to invoke more sophisticated constraints to observe 
different universal behavior from fc = geometric percolation [32| and so we expect the quantum behavior to be 
similar to quantum percolation for finite-dimensional lattices, though this conjecture should ultimately be tested. 

The paper is organized as follows. We review the methodology and results for quantum percolation on the cubic 
lattice [SOj as a means for calibration and then present our results for quantum fc = 0-core percolation and quantum 
fc = 3-core percolation on a Bethc-likc lattice. We conclude with a discussion of the implications of our results. 



II. QUANTUM PERCOLATION ON THE CUBIC LATTICE 



First, we analyze the level statistics for quantum percolation on the cubic lattice. To do so, we diagonalize the 
Hamiltonian defined by Eqs. 1 and 2 on the cubic lattice of length L with periodic boundary conditions to obtain a 



sequence of eigenvalues. This sequence is calculated for different realizations and increasing system sizes. The average 
density of states (DOS) for L = 15 as a function of occupation probabilities is presented in Fig. 2. The sharp peaks 
are due to small disconnected structures as discussed in Ref. [3^ . Now, one can apply the various measures of level 
statistics only if the density of states is smooth. There exists a smooth energy range around 0.4. The eigenenergies 
near this range are then arranged from highest to lowest and the nearest-neighbor level spacing, S, is calculated and 
subsequently normalized by the average nearest-neighbor level spacing, i.e. normalized level spacing is denoted by 
s = S/ <S>. 

One can then study the probability distribution for these level spacings, P{s), as a function of p and L. When 
the system is in the insulator regime, the eigenfunctions are localized and, therefore, do not interact with each other 
such that P{s) is Poisson distributed, i. e. P(s) = e"'*. When the system is in the metallic regime, Altshuler 
and Shklovskii [2g argued that if the width of the energy band of a sample E < Ec = hD/L^, where L'^/D is the 
characteristic time for an electron to diffuse through the sample, then the Hamiltonian of the system is characteristic 
of Gaussian Orthogonal Ensemble (GOE) in the absence of a magnetic field or spin-orbit scattering. More specifically, 
P{s) obeys the Wigner-Dyson distribution, P{s) = ^exp{—js'^). The cubic conductance < G >— j- < N{Ec) >= 

J- g tends to infinity when L — > cxd in the metallic regime, thus E < Ecis satisfied, and the level spacing distribution 
obeys the Wigner-Dyson distribution for the GOE. 

The plot of P{s) as function of bond occupation probability for L — 15 and an energy range of [0.2, 0.6] is displayed 
in Fig. 2. It can be seen that the expected transition from Wigner-like behavior for large p to Poisson behavior for 
small p. Note that all curves intersect at s ~ 2 as observed in the Anderson model \29\ . Other energy ranges studied 
yielded similar results. A convenient way to obtain the critical exponent v, characterizing the diverging localization 
length at the transition, is to study the parameter, 7, defined as 

which characterizes the transition from Wigner to Poisson as p is decreased. Note that 7 should increase from to 
1, as P{s) goes from Wigner to Poisson. Denoting ^(p) as the localization length (such that ^{p) ~ {p — Pq)""), 
the above parameter is expected to demonstrate scaling behavior, ^{p,L) — f[L/£,{p)]. In the vicinity of the critical 
quantum bond probability Pq , 

^{p,L)=^{pq,L) + C\^-l\L'/'', (4) 

Pq 

where C is a constant. 

Figure 3 plots 7 as a function of p around pq for increa sing L. The curves intersect near a single point given by 
Pg = 0.334. This result is consistent with the result in Ref. [30[. The single crossing point indicates that there exists a 
metal-insulator transition and one can apply the scaling collapse suggested above. The o ptim al scaling collapse shown 
in Fig. 3 yields v — 1.6 ±0.05, which should be compared to 1.32 ±0.08 obtained in Ref. [30|. Our latest result is even 



closer to the Anderson result than the previous work, where the most precise measurement is i^ = 1.58 ± 0.02 33| 



III. QUANTUM fc-CORE PERCOLATION ON A BETHE-LIKE LATTICE 

Following the procedure presented in Ref. [34], a Bethe-like lattice is numerically created by first considering a 
one-dimensional ring with N sites such that each lattice site has two bonds emanating from it. Next, additional 
bonds are constructed between different lattice sites at random. The number of random pairs connecting different 
lattice sites depends on the fixed coordination number, z. More specifically, there must be z bonds per site. See Fig. 
4 as an example for z = 3. As N increases, the average number of loops of length I increases as {z — 1)' such that 
the fraction of all lattice sites belonging to any loop of length < / for Z << \og{N) / \og{z — 1) is negligible. Therefore, 
the structure becomes increasingly tree-like as N increases and no surface effects that have to be dealt with, which 
occurs when one performs numerical simulations on the Bethe lattice. 

A. fc = 

For fc = 0, A. B. Harris J18l| gave a theoretical prediction of Pq on the Bethe lattice. One must simply solve 
1 + (pqcr^)^^ — (p<jCr)^^*-°'~^^ with a = z — 1. By solving Harris's self-consistency equation above for z = 6, Pq — 0.265. 
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FIG. 2: Left: The DOS for L — 15 and different bond occupation probabilities. Rigfit: P{s) for L = 15 for different ps. The 
Wigner and Poisson forms are also shown. 
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FIG. 3: Left: The function 'y{p,L) for different system sizes. Right: Scaling collapse for the cubic lattice. 



We can test this result numerically with the level statistics analysis. Figure 5 plots P{s) for different ps with N — 1024, 
z = 6 and an energy range of [0.3,0.8]. Figure 6 plots j{p,N) for z = 6. The curves intersect near p ~ 0.3. which 
is close to the analytical result for the Bethe lattice, though the agreement is not precise. This difference between 
the analytical calculation and our numerical calculation is presumably due to the nature of the lattice such that the 
quantum mechanics is much more sensitive to loops than geometric percolation. To test this notion, when measuring 
the onset of geometric percolation on the Bethe-like lattice, we do arrive at good agreement between the analytical 
result, Pc = l/(z — 1), and our numerical result. 

Moreover, the crossing point in Figure 6 indicates that we can collapse the data by assuming that, instead of 
j{p,L) — f[L/(^{p)] in the cubic lattice case, j{p,L) — f[N/N*{p)], where N* is a crossover size similar to the 
localization length in the three-dimensional case. In other words, N* ~ {p — pq)~^ . Therefore, j{p,L) = f[N^^'^ {p — 
Pq)]. It has been conjectured that v' = du^MF, where d^ is the upper critical dimension and vmf is the mean field 
correlation length exponent |;35|. For geometric percolation, d„ = 6 and i^mf — 1/2 such that v' = 3. To date, the 
upper critical dimension of quantum percolation is not known. Interestingly, the upper critical dimension for the 
Anderson model is potentially infinite such that the mean field correlation length scales with dimension d [3a] . More 
precisely, i' = 5 + ^• 

Figure 6 shows the optimal scaling collapse to yield the exponent i/' = 4.5 ± 0.2 assuming that Pq = 0.300(1). If we 
assume Pq = 0.265, then we do not arrive at a good scaling collapse. Regarding h'', if we assume the same mean field 
correlation length of vmf = 1/2 as in classical fc = percolation, which is consistent with the previous work of A. B. 
Harris ^18j showing that the mean field susceptibility for zero-energy eigenstates diverges with the same exponent as 
in the geometric percolation problem (with pq replacing pc), then we extract an upper critical dimension of (i„ = 9. 



If, on the other hand, the upper critical dimension is infinite as in the case of the Anderson model, then different 
analysis must be undertaken. 




FIG. 4: An example of a Bethe-like lattice with coordination number, z — 3. Dotted curves indicate the random pairs. 




FIG. 5: P{s) for the Bethedike lattice for N = 1024 and fc = for different ps. 



B. fc 



In Ref. ^], we gave an example of a MIT driven by fc = 3-core disorder on the Bethe lattice with z = 4. More 
precisely, the quantum conduction as a function of occupation probability p is a random first-order transition with 
Pq = Pc- To obtain this result, we assumed that the phase randomized between levels on the Bethe lattice, i.e. 
Id >> Ic- Here, we use level statistics to test the robustness of our prior results. Since the geometric critical 
percolation occupation probability pc is 8/9, which is close to 1, we choose z = 6, fc = 3, whose pc — 0.603 )Sx Y2fA l23| . 
Numerically, Figure 7 gives the fraction of remaining samples versus occupation probability, denoted by F , after 
applying the fc-core constraint recursively for z = 6, fc = 3. One can find that in a ideally infinite system size system, 
all occupied bonds are removed due to the fc-core constraint for p < 0.604, while p > 0.604, occupied bonds remain. 
Note that Pc= Pq = 0.604 is quite close to the analytical pc for z = 6 and fc = 3 for the classical physics. 

Figure 8 plots 7 as a function oip for different system sizes. As the system size increases, larger error bars are shown 
for p < Pc, indicating less remaining occupied samples due to the fc-core constraint. With no remaining occupied 
bonds for p < pc, the system is insulating in a trivial sense and not in the sense that 7—^1. What happens for 
p > Pc? The data indicates that 7 tends to as the system size increases for p > Pc, indicating that it is a metallic 
system. There is no crossing point above Pc such that pc must equal pq as we obtained previously. Note that we 
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FIG. 6: Left: The function y{p) for different system sizes. Right: The scalling collapse for 7(p) 



could have observed a crossing point somewhere above pc to indicate that Pq > Pc, but we do not. The data imply a 
discontinuous transition since 7 tends towards zero as N becomes large for p > Pq. 

IV. DISCUSSION 

We have applied level statistics for quantum percolation on three different systems. The results on cubic lattice 
are consistent with the results from Ref. [30]. For quantum percolation on a Bethe-like lattice, we find a threshold 
probability that differs somewhat from Harris's Bethe lattice result [la]. This difference is due to the sensitivity of 
loops in quantum mechanics occurring in the Bethe-like lattice. Moreover, a scaling collapse with system size N yields 
a new critical exponent i/' = 4.5(2). Assuming the correlation length exponent is vmf is 1/2, the same as the mean 
field result for geometric percolation, which is consistent with the Harris calculation [18| , the upper critical dimension 
du = 9- We have also demonstrated the robustness of our previous work for quantum k = 3-core percolation [4|. 
The level statistics here shows Pq = Pc, as before. In fact, for p < pc, all bonds are removed and the system is an 
insulator, while for p > pc, none of the bonds are removed and the system is conducting quantum mechanically. The 
data suggest that the system immediately goes to the Wigner-Dyson regime without going through a different regime 
at the transition, indicating a discontinuous transition, which also agrees with our previous work [4]. 

Therefore, our work provides an important counterexample for the Mott versus Anderson MIT paradigm, where 
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FIG. 7: Fraction remaining samples versus occupation probability after applying the fc-core constraint, F, on the Bethe-like 
lattice for z = 6, fc = 3. 
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FIG. 8: 7 as a function of p for different system sizes for the z = 6,k = 3 Bethe-like lattice. 



disorder-driven (Anderson) MITs are continuous and interaction-driven (Mott-Hubbard) MITs are discontinuous. 
Our counterexample is due to the correlations in the disorder as a result of the fc-core constraint. Correlations in the 
disorder have been shown previously to "complicate" matters. For instance, in one-dimensional wires with long-range, 
correlated disorder, there exists a MIT, which does not happen in the short-range, uncorrelated disorder case [37| . 

As for experimental implications, there are many experiments in the realm of classical transport on ordinary 
percolating systems. See, for example, Ref. [3a] for a recent one on nanowire composites. In addition, transport 
in undoped graphene is linked to classical electronic transport on percolating networks, though quantum effects are 
also relevant [39|. However, can a quantum k — 3-core percolation transition ever be realized? An experiment 
has already been conducted with a collection of silver quantum dots sitting atop of a Langmuir monolayer at room 
temperature [40]. As the interparticle spacing decreases by compressing the floating particles together, the electronic 
transport goes from hopping to tunnelling to ordinary metallic transport. The authors argue that disorder in the 
particle size and in the charging energy probably does not drive the transition and, instead, conjecture a possible 
first-order Mott transition at room temperature. However, in light of the analysis of the onset of classical conduction 
for k — 3-core, where k encodes the scalar aspect of local mechanical stability in particle packings [23[, we argue for 
a possible classical correlated percolation transition in conduction. Perhaps a quantum analog of this experiment can 
be realized in low-temperature packings of metallic nanoparticles to search for this new universality class of quantum 
disorder-driven MITs. 

JMS kindly acknowledges funding support from NSF-DMR-CAREER-0645373. 
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